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1 Introduction 



The theory of circle packings and, more generaUy, of circle patterns enjoys in recent years a fast de- 
velopment and a growing interest of specialists in complex analysis. The origin of this interest was 
connected with the Thurston's idea about approximating the Riemann mapping by circle packings, 
see [Tl|, [RS|. Since then the theory bifurcated to several subareas. One of them concentrates 
around the uniformization theorem of Koebe-Andreev-Thurston, and is dealing with circle pack- 
ing realizations of cell complexes of a prescribed combinatorics, rigidity properties, constructing 
hyperbolic 3-manifolds, etc |T|, |B§, Q. 

Another one is mainly dealing with approximation problems, and in this context it is advanta- 
geous to stick from the beginning with fixed regular combinatorics. The most popular are hexago- 
nal packings, for which the C°° convergence to the Riemann mapping was established by He and 
Schramm [HS[ . Similar results are available also for circle patterns with the combinatorics of the 
square grid introduced by Schramm 1^. It is also the context of regular patterns (more precisely, 
the two just mentioned classes thereof) where some progress was achieved in constructing discrete 
analogs of analytic functions (Doyle's spiralling hexagon packings [BDS| and their generalizations 
including the discrete analog of a quotient of Airy functions pH| ], discrete analogs of exp(2;) and 
erf (z) for the square grid circle patterns (S| , discrete versions of z" and log z for the same class of 
circle patterns |BP], [AB|). And it is again the context of regular patterns where the theory comes 
into interplay with the theory of integrable systems. Strictly speaking, only one instance of such 
an interplay is well-established up to now: namely, Schramm's equation describing the square grid 
circle packings in terms of Mobius invariants turns out to coincide with the stationary Hirota's 
equation, known to be integrable, see |BP], [^. It should be said that, generally, the subject of 
discrete integrable systems on lattices different from Z" is underdeveloped at present. The list of 
relevant pubhcations is almost exhausted by ^D|, fNS], @, [IOPI . 

The present paper contributes to several of the above mentioned issues: we introduce a new 
interesting class of circle patterns, and relate them to integrable systems. Besides, for this class we 



construct, in parallel to [BP], [|AB| , the analogs of the analytic functions z", logz. 

This class is constituted by hexagonal circle patterns, or, in other words, by circle patterns with 
the combinatorics of the regular hexagonal lattice (the honeycomb lattice). This means that each 
elementary hexagon of the honeycomb lattice corresponds to a circle, and each common vertex 
of two hexagons corresponds to an intersection point of the corresponding circles. In particular, 
each circle carries six intersection points with six neighboring circles. Since at each vertex of the 
honeycomb lattice there meet three elementary hexagons, there follows that at each intersection 
point there meet three circles. 

This class of hexagonal circle patterns is still too wide to be manageable, but it includes several 
very interesting subclasses, leading to integrable systems. For example, one can prescribe inter- 
section angles of the circles. This situation will be considered in a subsequent publication. In the 
present one we consider the following requirement: the six intersection points on each circle have 
the multi-ratio equal to —1, where the multi-ratio is a natural generalization of the notion of a 
cross-ratio of four points on a plane. 

We show that, adding to the intersection points of the circles their centers, one embeds hexagonal 
circle patterns with the multi-ratio property into an integrable system on the regular triangular 
lattice. Each solution of this latter system describes a peculiar geometrical construction: it consists 
of three triangulations of the plane, such that the corresponding elementary triangles in all three 
tilings are similar. Moreover, given one such tiling, one can reconstruct the other two almost 
uniquely (up to an affine transformation). If one of the tilings comes from the hexagonal circle 
pattern, so do the other two. This results are contained in Sect. ^ ^. In the intermediate Sect. |3| 
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we discuss a general notion of integrable systems on graphs as flat connections with the values in 
loop groups. It should be noticed that closely related integrable equations (albeit on the standard 
grid I?) were previously introduced by Nijhoff R in a totally different context (discrete Bussinesq 



equation), see also similar results in [BK|. However, these results did not go beyond writing down 
the equations: geometrical structures behind the equations were not discussed in these papers. 

Having included hexagonal circle patterns with the multi-ratio property into the framework 
of the theory of integrable systems, we get an opportunity of applying the immense machinery 
of the latter to studying the properties of the former. This is illustrated in Sect. |5|, ^ where 
we introduce and study some isomonodromic solutions of our integrable system on the triangular 
lattice, as well as the corresponding circle patterns. Finally, in Sect. we define a subclass of these 
"isomonodromic circle patterns" which are natural discrete versions of the analytic functions z", 
logz. The results of Sect. constitute an extension to the present, somewhat more intricate, 
situation of the similar constructions for Schramm's circle patterns with the combinatorics of the 
square grid [AB|. 



2 Hexagonal circle patterns 



Figure 1: The regular triangular lattice with its hexagonal sublattices. 



First of all we define the regular triangular lattice TC as the cell complex whose vertices are 

V{TC) = ^i = k + £u; + mLo'^ : k,£,mezY where a; = exp(27ri/3), (2.1) 

whose edges are all non-ordered pairs 

E{TC) = {[31,32] : 31,32 e V{TC), ^ -32I = l}, (2.2) 

and whose 2-cells are all regular triangles with the vertices in V{TC) and the edges in E(TC). 
We shall use triples {k, i, m) G I? as coordinates of the vertices of the regular triangular lattice, 
identifying two such triples iff they differ by the vector (n, n, n) with n £ Z. We call two points 
31,32 neighbors in TC, iff [31,32] G E{T£). 

To the complex TC there correspond three regular hexagonal sublattices Ti-Cj, j = 0, 1,2. 
Each TCCj is the cell complex whose vertices are 

V{nCj) = \i = k + eu; + mu;'^ : k,e,meZ, k + e + m^j (mod 3)|, (2.3) 
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whose edges are 

E{nCj) = {[31,32] : 31,32 G V{nCj), |3i -32I = l}, (2.4) 

and whose 2-cells are all regular hexagons with the vertices in V{'HCj) and the edges in E{HCj). 
Again, we call two points 31,32 neighbors in TiCj, iff [31,32] G E(HCj). Obviously, every point 
in V{HCj) has three neighbors in TCCj, as well as three neighbors in TC which do not belong to 
V{HCj). The centers of 2-cells of TCCj are exactly the points of V{TL) \ V{J-LCj), i.e. the points 
l' = k + iuj + mu;^ with k + i + m = j (mod 3) . 

In the following definition we consider only HCq, since, clearly, HC\ and HC2 are obtained 
from HCq via shifting all the corresponding objects by w, resp. by 00^. 

Definition 1 We say that a map w : V{HjCo) <—>■ C defines a hexagonal circle pattern, if the 
following condition is satisfied: 

• Let 

= + £^ ^ V{HCo), fc = l,2, ...,6, where e = exp(7ri/3), 

be the vertices of any elementary hexagon in HCq with the center G V{TC)\V{7iCo). Then 
the points w{^i), 10(22), ■■ ■ ywide) G C lie on a circle, and their circular order is just the listed 
one. We denote the circle through the points u;(3i), 10(32), . • • ,'f«(36) by C{j,'), thus putting it 
into a correspondence with the center 3' of the elementary hexagon above. 

As a consequence of this condition, we see that if two elementary hexagons of HCq with the centers 
in l' il" G V{TC) \ V{TLCq) have a common edge [31,32] G E{Ti.Co), then the circles C(3') and C(3") 
intersect in the points ^(31), w{$2)- Similarly, if three elementary hexagons of HCq with the centers 
in 3', 3", 3'" G V{T C) \ V{nCo) meet in one point 30 G ViHCo), then the circles C{}'), C{i") and 
C{i"') also have a common intersection point u;(3o)- (Note that in every point 30 G V{J-LCq) there 
meet three distinct elementary hexagons of HjCq)- 

Remark. Sometimes it will be convenient to consider circle patterns defined not on the whole 
of HCq, but rather on some connected subgraph of the regular hexagonal lattice. 

We shall study in this paper a subclass of hexagonal circle patterns satisfying an additional 
condition. We need the following generalization of the notion of cross-ratio. 

Definition 2 Given a {2p) -tuple {wi,W2,... ,W2p) G C^^ of complex numbers, their multi-ratio 
is the following number: 

^ nj=i(^^2j-i -u'2j) 
M{wi,W2,... ,W2p) = — 7 7> (2-5) 

where it is agreed that u'2p+i = wi . 
In particular, 

X jwi - W2)iw3 - W4) 
M{Wi,W2,Ws,W4) = 

{W2 - W3){W4 - Wi) 

is the usual cross-ratio, while in the present paper we shall be mainly dealing with 

{Wl - W2){W3 - Wa){w^ - Wq) 
M{wi,W2,... ,W6) = -. 

[W2 - W3){W4, - W5)[W6 - Wi) 

The following two obvious properties of the multi-ratio will be important for us: 
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(i) The multi-ratio M{wi,W2,--- ,W2p) is invariant with respect to the action of an arbitrary 
Mobius transformation w i— > {aw + b)/{cw + d) on all of its arguments. 

(ii) The multi-ratio M{wi,W2, ■ ■ ■ ,W2p) is a Mobius transformation with respect to each one of 
its arguments. 

We shall need also the following, slightly less obvious, property: 

(iii) If the points 101,1112, ■■ ■ , W2p-i lie on a circle C C C, and the multi-ratio M{wi, W2, ■ ■ ■ , W2p) 
is real, then also W2p G C. 



Definition 3 We say that a map w : V{HCo) i—^ C defines a hexagonal circle pattern with 

MR = — 1, if in addition to the condition of Definition [J the following one is satisfied: 

• For any elementary hexagon in HCq with the vertices 31,325 •• • )36 S V{7iCo) (listed counter- 
clockwise), the multi-ratio 

M{wi,W2,... ,we) = -l, (2.6) 

where Wk = u^(3fc)- 



Geometrically the condition ( |2.6D means that, first, the lengths of the sides of the hexagon with 
the vertices W1W2 ■ ■ -Wq satisfy the condition 

\wi — W2\ ■ \W3 — W^l ■ \W5 — Wq\ = \W2 — W^l ■ \W4 — W^l ■ \wq — Wi\, 

and, second, that the sum of the angles of the hexagon at the vertices wi, W3, and W5 is equal 
to 27r(mod 27r), as well as the sum of the angles at the vertices W2, W4, and wq. Notice that if 
a hexagon is inscribed in a circle and satisfies (|2.6| ), then it is conformally symmetric, i.e. there 
exists a Mobius transformation mapping it onto a centrally symmetric hexagon. Notice also that 
the regular hexagons satisfy this condition. 

To demonstrate quickly the existence of hexagonal circle patterns with MR = — 1 we give their 
construction via solving a suitable Cauchy problem. 

Lemma 4 Consider a row of elementary hexagons of TiCo running from the north-west to the 
south-east, with the centers in the points 3^ = k — koo. Let the map w be defined in five vertices 
of each hexagon ~ in all except 3'^ -|- e. Suppose that the five points 1^(3'^ + e-'), j = 2, 3, . . . , 6, lie 
on the circles C(3^). These data determine uniquely a map w : V{T-LCq) ^ C yielding a hexagonal 
circle pattern with MR = —1 on the whole lattice. 



Proof. Equation (2^) determines the points ti'(3^ -|- s), which, according to the property above, 
lie also on C(3^). Now for every hexagon of the parallel row next to north-east, with the centers in 
the points 3'^' = ^'f^ + 1 + e = {k -\- 2) — {k — l)io, we know the value of the map vj in three vertices, 
namely in 

3fc + = 3fc + 1 = 3fc+i + dk + = d'k + ^; ik + = d'k+i + 

This uniquely defines the circle C {],'[), as the only circle through three points ^(3'^' -|- e^), w{i'l. + e^) 
and w{i'l. -\- e^). The intersection points of these circles of the second row give us the values of the 
map w in the points 3'^' -|- and 3'^' -|- . Namely, w{i'l + e^) is the intersection point of C(3'^') with 
C(3fc_i)) different from ^(3'^' + e^), and w{i'l^ + e^) is the intersection point of C(3'^') with C{i'l_^-^), 
different from w{i'l-\- £^). Therefore we get the values of the map w in five vertices of each hexagon 



5 



of the next parallel row - in all except 3^ + e. The induction allows to continue the construction 
ad infinitum. I 

Now we show that, adding the centers of the circles of a hexagonal pattern with MR = — 1 to 
their intersection points, we come to a new interesting notion. 

Theorem 5 Let the map w : V{7iCQ) 1— > C define a hexagonal circle pattern with MR = — 1. 
Extend w to the points ofV{TC) \ V{7iCQ) by the following rule. Fix some point Pqo £ C. Let 3' be 
a center of an elementary hexagon ofTiC^. Setw{])') to be the reflection of the point in the circle 
C{i'). Then the condition ( ^.(\ ) holds also for = w{}k) the case when the points 31,32; • • • ,36 
are the vertices of any elementary hexagon of the two complementary hexagonal sublattices TiCi 
and HC2. 




Figure 2: An elementary hexagon with its center point sent to 00. 



Proof. Consider the situation corresponding to an elementary hexagon of the sublattice TiCi or 
Ti.C2 (see Fig. ^). The point wq is the intersection point of the three circles C(3i), 0(33), and 
(7(35), the points wi, W3, and are obtained by reflection of Poo in the corresponding circles, and 
the points W2, W4, and wq are the pairwise intersection points of these circles different from wq. 
To simplify the geometry behind this situation, perform a Mobius transformation sending wq to 
infinity. Then the circles C(3i), (7(33), and (7(35) become straight lines, and the points wi, W3, 
are the reflections of Poo in these lines (see Fig. ^; for definiteness we suppose here that the Mobius 
image of Poo lies in the interior of the triangle formed by these straight lines). By construction, 
one gets: 

\W2 — Wi\ = \W2 — Ifsl, |l«4 — = |tU4 — W^\, \w% — = \wq — Wi\\ 

the angles by the vertices W2, w^, wq are equal to 2(ai + 02), 2(/?i + P2), 2(71 + 72), respectively, 
so that their sum is equal to 

2(ai + a2 + Pi+ P2 + 71 + 72) = 27r; 



6 



the angles by the vertices wi, W3, are equal to vr — (ai + 72), vr — + 02), vr — (71 + 
respectively, so that their sum is equal to 

Svr - (ai + a2 + Pi + P2 + li + 72) = 27r. 
This proves that the hexagon under consideration satisfies ( p.GI ). I 

A particular case of the construction of Theorem |5| is when Pqo = 00, so that the map w is 
extended by the centers of the corresponding circles. In any case, this theorem suggests to consider 
the class of maps described in the following definition. 

Definition 6 We say that the map w : V{TC) C defines a triangular lattice with MR = 

— 1, if the equation ( \2.(\ ) holds for Wk = w{^k), whenever the points 31,321 • • • ,36 o'^e the vertices 
(listed counterclockwise) of any elementary hexagon of any of the sublattices Ti-Cj {j = 0, 1, 2). 

In the next section we shall discuss an integrable system on the regular triangular lattice, each 
solution of which delivers, in a single construction, three different triangular lattices with MR = — 1. 
However, these three lattices are not independent: given such a lattice, the two associated ones 
can be constructed almost uniquely (up to an affine transformation w aw + h). It will turn out 
that if the original lattice comes from a hexagonal circle pattern with MR = —1, then the two 
associated ones do likewise. 

3 Discrete flat connections on graphs 

Let us describe a general construction of "integrable systems" on graphs which does not hang on 
the specific features of the regular triangular lattice. This notion includes the following ingredients: 

• An oriented graph G; the set of its vertices will be denoted V{Q), the set of its edges will be 
denoted E{Q). 

• A loop group G[A], whose elements are functions from C into some group G. The complex 
argument A of these functions is known in the theory of integrable systems as the "spectral 
parameter" . 

• A "wave function" ^ : V{Q) ^ ^[A], defined on the vertices of Q. 

• A collection of "transition matrices" L : E{Q) ^ G[A] defined on the edges of Q. 

It is supposed that for any oriented edge e = (31,32) G E{Q) the values of the wave functions in its 
ends are connected via 

^'(32,A)=L(e,A)*(3i,A). (3.1) 

Therefore the following discrete zero curvature condition is supposed to be satisfied. Consider any 
closed contour consisting of a finite number of edges of Q: 

ei = (31)32), 22 = (32,33), ep = (3p,3l)- 

Then 

L(ep,A)---L(e2,A)L(ei,A) = /. (3.2) 
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In particular, for any edge e = (31,32)) if e ^ = (32,31), then 

L(e-\A) = (L(e,A))"\ (3.3) 

Actually, in applications the matrices L{c, A) depend also on a point of some set X (the "phase 
space" of an integrable system), so that some elements x(e) G X are attached to the edges e of 
G- In this case the discrete zero curvature condition (|3.2| ) becomes equivalent to the collection of 
equations relating the fields x(ei), . . . , x(cp) attached to the edges of each closed contour. We say 
that this collection of equations admits a zero curvature representation. 

For an arbitrary graph, the analytical consequences of the zero curvature representation for 
a given collection of equations are not clear. However, in case of regular lattices, like TC, such 
representation may be used to determine conserved quantities for suitably defined Cauchy problems, 
as well as to apply powerful analytical methods for finding concrete solutions. 

Remark. The above construction of integrable systems on graphs is not the only possible one. 
For example, in the construction by Adler @ the fields are defined on the vertices of a planar 
graph, and the equations relate the fields on stars consisting of the edges incident to each single 
vertex, rather than the fields on closed contours. Examples are given by discrete time systems of 
the relativistic Toda type. In the corresponding zero curvature representation the wave functions ^ 
naturally live on 2-cells rather than on vertices. The transition matrices live on edges: the matrix 
L(e, A) corresponds to the transition across e and depends on the fields sitting on two ends of e. 

4 An integrable system on the regular triangular lattice 

We now introduce an orientation of the edges of the regular triangular lattice TC. Namely, we 
declare as positively oriented all edges of the types 

(3,3 + 1), (3,3 +'^), (3,3 +'^^)- 
Correspondingly, all edges of the types 

(3,3-1), (3,3 -cj), (3,3 -w^) 

are negatively oriented. Thus all elementary triangles become oriented. There are two types of 
elementary triangles: those "pointing upwards" (3,3 + 1^,3 — 1) are oriented counterclockwise, while 
those "pointing downwards" (3,3 +cl'^,3 — 1) are oriented clockwise. 

4.1 Lax representation 

The group G[X] we use in our construction is the twisted loop group over SL(3,C): 



L:Ch^SL(3,C) L(a;A) = J1L(A)1^"^|, (4.1) 



where Q = diag(l, cj, a;^). The elements of G[X] we attach to every positively oriented edge of TC 
are of the form 

1 A/ 

L(A) = (1 + A3)-i/3 I 1 Xg\, fgh=l. (4.2) 

Xh 1 
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Hence, to each positively oriented edge we assign a triple of complex numbers (/, h) G satis- 
fying an additional condition fgh = 1. In other words, choosing {f,g) (say) as the basic variables, 
we can assume that the "phase space" X mentioned in the previous section, is C* x C,,. The scalar 
factor (1 + A'^)"^/'^ is not very essential and assures merely that detL(A) = 1. 



It is obvious that the zero curvature condition (3^) is fulfilled for every closed contour in TC, 
if and only if it holds for all elementary triangles. 

Theorem 7 Let ei, e2, es be the consecutive positively oriented edges of an elementary triangle of 
TC. Then the zero curvature condition 

L(e3,A)L(e2,A)L(ei,A) =/ 

is equivalent to the following set of equations: 

/i + /2 + /3 = 0, 51 + 52 + 53 = 0, (4.3) 

and 

fi9i = /352 /252 = /i53 <^ /353 = /25i, (4-4) 

with the understanding that hj. = {fkOk) "^) k = 1,2,3. 

Proof. An easy calculation shows that the matrix equation L3L2L1 = I consists of the following 
nine scalar equations: 

/i + /2 + /3 = 0, 51 + 52 + 53 = 0, /ii + /i2 + /13 = 0, (4.5) 
/352/ii = 1, 53^2/1 = 1, /i3/25i = 1, (4.6) 

/352 + /351 + /251 = 0, 53/12 + 53/ii + 52/11 = 0, /13/2 + hsfi + /i2/i = 0. (4.7) 

It remains to isolate the independent ones among these nine equations. First of all, equations ( [4.7| ) 
are equivalent to (^), provided (^) and fudk^k = 1 hold. For example: 

73(52 + 51) + /251 = <^ /353 = /251 <^ /13/251 = 1- 
Next, the conditions fkQkhk = 1 allow us to rewrite ( |4.6D as 

/i5i = /352, /252 = /l53, /353 = /25l- (4.8) 
Further, all equations in (^^) are equivalent provided ( |4.3| ) holds. For example: 

/i5i = /352 ^ (/2 + /3)5i = 73(51 + 53) =^ /251 = /353- 
Finally, hi + h2 + h^ = ^ follows from (U), (gj). Indeed, 

/11 + /12 = (/i5i)~' + (/252)-' = (/352)-' + (/252)-' 

= (/252)-'(/2 + /3)/3"' = -(/252)-Vl/3'' 
= -(/l53)^Vl/3"' = -(/353)-' = -/13. 

The theorem is proved. For want of a better name we shall call the system of equations ( [1.3D , ( [4.4| ) 
the fgh system. I 



9 



The equations ( [4. 5]) may be interpreted in the following way: there exist functions u,v,w : 
V{TC) I— > C such that for any positively oriented edge e = (31,32) there holds: 

/(e) = n(32) - u(3i), g{c) = v{i2) - v{u), H'^) = w{i2) - w{u)- (4.9) 

The function u is determined by / uniquely, up to an additive constant, and similarly for the 
functions v, w. Having introduced functions u, v, w sitting in the vertices of TC, we may reformulate 
the remaining equations ( [4.4| ) as follows: let 31,32,33 be the consecutive vertices of a positively 
oriented elementary triangle, then 

^(32) -^(31) ^ ^'(33) -^(32) 

The equations arising by cyclic permutations of indices (1,2,3) ^ (2,3,1) are equivalent to this 
one due to ( |4.4| ). So, we have one equation pro elementary triangle 313233- Its geometrical meaning 
is the following: the triangle u{},i)u{}2)u{i3) is similar to the triangle f (32)^(33)^(31) (where the 
corresponding vertices are listed on the corresponding places). Of course, these two triangles are 
also similar to the third one, w{-^3)w{},i)w{-^2)- 

4.2 Cauchy problem 

We discuss now the Cauchy data which allow one to determine a solution of the fgh-system. The 
key observation is the following. 

Lemma 8 Given the values of two fields, say u and v, in three points 30, 31 = 30 + I CLnd}2 = 30+*^; 
the equations of the fgh-system determine uniquely the values ofu andv in the point ^3 = 30 + l+aj; 

U3 - uo = [ui - uo) \-{u2-uo) , (4.11) 

Vi -V2 V2 - Vi 

^^3 - ■'^i = (?^i - ■Wo) ^ -f^s - ■"2 = (^^2 - I'd) • (4.12) 

Uo - U3 Uq- U3 



Proof. The formula ( 4.11 ) follows by eliminating V3 from 

Up - U3 _ Vl - Vq U0-U3 _ V2- vq 

Ui - Uq V3-V1'' U2- Uq V3 - V2 

These equations yield then ( 4.12| ). I 



(4.13) 



This immediately yields the following statement. 

Proposition 9 a) The values of the fields u and v in the vertices of the zig-zag line running 
from the north-west to the south-east, 

^} = k + euj:k + i = 0,l|, 

uniquely determine the functions u, v : V{TC) ^ <C on the whole lattice. 

b) The values of the fields u and v on the two positive semi-axes, 

= /c : A; > o| U |3 = ^cj : £ > o|, 

uniquely determine the functions u, v on the whole sector 

= k + iuj : k,i > 0^ = e V{TC) : < arg(3) < 27r/3| . 

Proof follows by induction with the help of the formulas ( [4.11 ), ( 4.12| ). I 
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4.3 Sym formula and related results 



There holds the following result having many analogs in the differential geometry described by 
integrable systems ("Sym formula", see, e.g., [[BP| ) . 



Proposition 10 Let ^(3, A) be the solution of ( \3. 1\ ) with the initial condition ^'(jo. A) = / for some 
30 G V{TC). Then the fields u,v,w may he found as 



d^ 





( 


u 
















\ w 





/ 



(4.14) 



Proof. Note, first of ah, that from ^'(jo, 0) = / and L(e, 0) = / there follows that ^(3, 0) = I for 
all 3 € V{TC). Consider an arbitrary positively oriented edge e = (31,32)- From (^) there follows: 



d^in) d^in) 



dX 



dX 



'^(^)^(30+i(e)'''^^^^^ 



dX 



dX 



At A = we find: 



^^(32) 




dL[t) 




dX 


A=0 dX 


A=0 


A=0 ( 




' 


An) - u{ix) 










^^(32) 




V ^^(32) - w{j>\) 








dX 



/(e) 
g{z) 
h{z) 



- ^(31^ 




This proves the Proposition. I 

Next terms of the power series expansion of the wave function ^'(3, A) around A = also deliver 
interesting and important results. 



Proposition 11 Let ^'(3, A) be the solution of (3.1) with the initial condition ^'(39, A) = / for some 
30 e V{TC). Then 



1 














b 









\ 


c 


/ 



2 dX^ 

where the function a : V{TC) ^ C satisfies the difference equation 

0(32) - a(3i) = ^^(31) (-"(32) - ^^(31)) , 



(4.15) 



(4.16) 



and similar equations hold for the functions b,c : V{TC) ^ C (with the cyclic permutation 
{u,v,w) I— >• (w,u,v)). 



Proof. Proceeding as in the proof of Proposition we have: 



d^v&(32) ^^^(31) _ (d^Ljc) 
dA2 dA2 V dX"^ 



-^{n) + 2 



dLjc) d^in) 
dX dX 



+ L(e)- 



dP-^{ll)\ d2^(3i) 



dX^ 



dX^ 
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Taking into account that <f'L{z)/d}?\\=Q 



0, we find at A = 0: 




This imphes the statement of the proposition. I 

Notice that it is a priori not obvious that the equation ( 4.16| ) admits a well-defined solution on 
V{TC), or, in other words, that its right-hand side defines a closed form on TC. This fact might 
be proved by a direct calculation, based upon the equations of the f gh-system, but the above 
argument gives a more conceptual and a much shorter proof. 



Corollary 12 Under the conditions of Propositions we have: 

1 



+ 



A=o 



d"^ 



2 

A=0 




2 dA2 

where the function a' : V(TC) C satisfies the difference equation 

a' (32) - a' (31) = ^^(32) (^(32) - 1^(31)) , 



(4.17) 



and similar equations hold for the functions b',c' : V(TC) 
(m, V, w) 1-^ {w, u, v) ). 



(4.18) 

C (with the cyclic permutation 



Further examples of such exact forms may be obtained from the values of higher derivatives of the 
wave function ^'(3, A) at A = 0. 



4.4 One— field equations 

We discuss now the equations satisfied by the field u alone, as well as by the field v alone. In this 
point we make contact with the geometric considerations of Sect. |2[ 

Theorem 13 1. Both maps u,v : V(TC) 1— > C define triangular lattices with MR = — 1. 
In other words, if u^i2,--- ,}6 cl^^ the vertices (listed counterclockwise) of any elementary 
hexagon of any of the hexagonal sublattices Ti-Cj (j =0,1,2), and ifuj. = ^(3^) and Vf^ = ^(3^), 
then there hold both the equations 



M{ui,U2, ... ,ue) 



(4.19) 



and 



M{vi,V2, ... ,^6) 



-1. 



(4.20) 



2. Given a triangular lattice u : V{TC) ^ C with MR = —\, there exists a unique, up to an 
affine transformation v ^ av + b, function v : V{TC) 1— > C such that are satisfied 

everywhere. This function also defines a triangular lattice with MR = —1. 
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3. Given a pair of complex-valued functions {u, v) defined on V{TC) and satisfying the equa- 
tion ^-lOj ) everywhere, there exists a unique, up to an affine transformation, function w : 
V{TC) I— > C such that the pairs {v,w) and {w,u) satisfy the same equation. The function w 
also defines a triangular lattice with MR = — 1. 

Proof. 1. To prove the first statement, we proceed as follows. Let 3' £ V{TC), and let the vertices 
of an elementary hexagonal with the center in 3' be enumerated as 3^ = 3' + s^, k = 1,2, .. . ,6. 
Then the following elementary triangles are positively oriented: (32^532^-1 ,3') and {12k ■,l2k+i ■,]>') for 
= 1,2,3 (with the agreement that 37 = 31). According to ( [4.10[ ), we have: 

U2k-l - U2k _ V' - V2k-1 U2k+1 - U2k _ v' - V2k+1 r, _ r, o 

U' - U2k-1 V2k -V' U' - U2k+1 V2k " V 

Dividing the first equation by the second one and taking the product over A; = 1, 2, 3, we find: 

3 

TT- U2k-l — U2k _ ^ 
^J-^ U2k+1 - U2k 



which is nothing but ( 4.19 ). The proof of ( [1.20D is similar. 



2. As for the second statement, suppose we are given a function u on the whole of V{T£,). For 
an arbitrary elementary triangle, if the values of v in two vertices are known, the equation ( [4.4| ) 
allows us to calculate the value of v in the third vertex. Therefore, choosing arbitrarily the values 
of V in two neighboring vertices, we can extend this function on the whole of V{TC), provided 
this procedure is consistent. It is easy to understand that it is enough to verify the consistency in 
running once around a vertex. But this is assured exactly by the equation ( 4.19| ). 



3. To prove the third statement, notice that the proof of Theorem |^ shows that the formula 

1 1 

h{e) = w{i2) - w{ii) = ——-- = —— , . . (4.21) 

/(e)c/(e) (u(32) - u{u)){v{i2) - v{u)) 

valid for every edge e = (31,32) of TC, correctly defines the third field h of the /g/i-system. All 
affine transformations of the field w thus obtained, and only they, lead to pairs (v, w) and (w, u) 
satisfying (|1^). ■ 

Remark. Notice that the above results remain valid in the more general context, when the 
fields /, g, h do not commute anymore, e.g. when they take values in H, the field of quaternions. 
The formulation and the proof of Theorem hold in this case literally, while the formula ( |4.19| ) 
reads then as 

{Ui - U2){U2 - U^Y^{u^ - Ui){Ui - U5)~^{U5 - U6){uq - Ui)~^ = -1, (4.22) 

and similarly for v,w. 
4.5 Circularity 

Recall that hexagonal circle patterns with MR = — 1 lead to a subclass of triangular lattices with 
MR = — 1, namely those where the points of one of the three hexagonal sublattices lie on circles. 
We now prove a remarkable statement, assuring that this subclass is stable with respect to the 



transformation u ^ v described in Theorem 13 
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Theorem 14 Letu : V{TCCj) C define a hexagonal circle pattern with MR = —1. Extend it with 
the centers of the circles to u : V['TC) i— > C, a triangular lattice with MR = —1. Letv : V{TC) i— > C 



be the triangular lattice with MR = —1 related to u via ( ^.IC ). Then the restriction of the map v 
to the sublattice TCCj^i also defines a hexagonal circle pattern with MR = —1, while the points v 
corresponding to TC \ TiCjj^i are the centers of the corresponding circles. 



Proof starts as the proof of Theorem 13. Let 3' be a center of an arbitrary elementary hexagon of 



the sublattice HCj^i, i.e. 3' = k + iio + mu)'^ with k + £ + m = j + 1 (mod 3). Denote by 3^ = 3' + e'^, 
fe = 1, 2, . . . ,6 the vertices of the hexagon. As before, considering the positively oriented triangles 
{d2k,d2k-i,d') and (32A;,32fe+i,3')! k = 1,2,3, surrounding the point 3', we come to the relations 

U2k-1-U2k v'-V2k-l U2k+1-U2k v' - V2k+1 , . „„ . , 

— = -, — = fe = l,2,3. (4.23) 

U' - U2k~l V2k -V' v! - U2k+1 V2k " V 

But, obviously, 32^-1 {k = 1,2,3) are centers of elementary hexagons of the sublattice 'HCj. By 
condition, the points U2k-2i U2k and u' lie on a circle with the center in U2k-i- Therefore, 

\u2k - U2k-i\ = \u2k-2 - U2k-i\ = \u' - U2k~i\, k = 1,2,3. (4.24) 



So, the absolute values of the left-hand sides of all equations in ( 4.23| ) are equal to 1. It follows 
that all six points wi, f2, . . . , lie on a circle with the center in v' . I 



5 Isomonodromic solutions 

Recall that we use triples (fc, £, m) £ 1? as coordinates of the vertices 3 = /c + ^u; + muP' , and that 
two such triples are identified iff they differ by the vector (n, n, n) with n G Z. By the /c-axis we 
call the straight line R C C, resp. by the -^-axis the straight line Ruj, and by the m-axis the straight 
line Ruj'^. 

It will be sometimes convenient to use the symbols ~, " and ~ to denote the shifts of various 
objects in the positive direction of the axes k, i, m, respectively, and the symbols ^, ^, to denote 
the shifts in the negative directions. This will apply to vertices, edges and elementary triangles of 
TC, as well as to various objects assigned to them. For example, if 3 G V{TC), then 

I=3 + l> ^ = 3-1, ?=3+'^> i=}-^^, d=2 + ^'^, 3=3-0;^. 

Similarly, if e = (31,32) E E(TC), then 

1= (31 + 1,32 + 1), ?= (31 +a;,32 +t^), e = (31 + w^, 32 + w^), etc. 

A fundamental role in the subsequent presentation will be played by a non- autonomous con- 
straint for the solutions of the fgh-system. This constraint consists of a pair of equations which are 
formulated for every vertex 3 G V{TC) and include the values of the fields on the edges incident to 
3, i.e. on the star of this vertex. It will be convenient to fix a numeration of these edges as follows: 

eo = (3,3), e2 = (3,?), e4 = (3,3), (5.1) 
ei = (3,3), e3 = (l,3), 25 = (^,3)- (5-2) 

The notations /o, • • • , /e will refer to the values of the field / on these edges: 

/o = u - n, f2 = u-u, f^ = u-u, (5.3) 
fi = u- u, f3 = u- u, f^ = u- u, (5.4) 
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Figure 3: Notations for u and /. 



and similarly for the fields g, h, see Fig. ^. 
The constraint looks as follows: 

_ ^ /ogo/3 72^2/5 |_ ^ f 494/1 

fo9o + 90/3 + f393 f292 + 92/5 + /sffs /454 + g^fi + fi9i ' 

p^. - J, aof-m ^ ^ 92h95 ^ gdm ^^ g^ 

fo90 + 50/3 + 7353 /252 + 52/5 + /555 /454 + 94fl + fl9l' 

These are supposed to be the equations for the vertex 3 = /c + £cj + mw^, and we use the notations 
u = u{}), V = v{}). Since the fields u, v are defined only up to an affine transformation, one should 
replace the left-hand sides of the above equations by au + 0, f3v + ip, respectively, with arbitrary 
constants (j), ip. In the form we have choosen it is imposed that the fields u, v are normalized to 
vanish in the origin. 

Proposition 15 The equations (\5.3^), ( |5. are well defined equations for the point 3 G V(TC), 
i.e. they are invariant under the shift {k,£,m) 1— > (k + n,£ + n,m + n), provided the equations i\4-l(\) 
hold. 

Proof is technical and is given in the Appendix I 

We mention an important consequence of this proposition. Apparently, the constraint (^.5|), 
(|5.6| ) relates the values of the fields u, v in seven points shown on Fig. |3[ However, we are free to 
choose any representative {k, i, m) for 3. In particular, we can let vanish any one of the coordinates 
k, i, m. In the corresponding representation the constraint relates the values of the fields u, v in 
five points, belonging to any one of the three possible four-leg crosses through 3. 

An essential algebraic property of the constraint (|5.5| ), ( |5.6D is given by the following statement. 



Proposition 16 // the equations hold, then the constraints ( |5.5[ j, ( ^5. 6|j imply a similar 

equation for the field w (vanishing at i = {)): 

^ fogo + 50/3 + /353 ^ f292 + 52/5 + f595 /454 + 54/1 + fl9l ' 

where 7 = 1 — a — 

Proof is again based on calculations and is relegated to the Appendix p|. I 
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Remark. We notice that restoring the fields hk = ^/{fkdk) allows us to rewrite the equations 
H), (13) as 

PV = k ; ; h« ; ; hm ; ; — , (5.8) 

90^0 + hogs + gshs 52/12 + "-255 + 55% 94/14 + /i45'i + gihi 

= k I £ i m ^'^"^^ (5.9) 

hofo + fohs + hsfs /12/2 + /2/i5 + /14/4 + fihi + hifi ' 



which coincides with ( |5.5| ) via a cyclic permutation of fields (/, g, h) ^ (g, h, f) performed once or 
twice, respectively, and accompanied by changing a to /?, 7, respectively. 

Another similar remark: as it follows from the formulas (B.2), (|B.3|) used in the proof of 



Proposition |l5| (and their analogs for the fields g, h), the constraints (|5.5| ), (|5.6| ), ( |5.7D may be 
rewritten as equations for the single field n, resp. v, w: 

_ , fofsifi + 12) g 72/5 (/a + /4) fifiifb + /o) (^^(^'\ 

" (/0-/2)(/l-/3)+ (/2-/4)(/3-/5)+'"(/4-/o)(/5-/l)' ^ ^ 

^ ^„ g0g3(gl+g2) ff2g5(g3+g4) 3491(35+ go) 

(go -g2)(gl -gs) (g2 - g4)(g3 - g5) (g4 -go)(g5 -gl)' 

_ hohs{hi + h2) /i2%(% + %) h^hiQi^ + /ip) 

" (/10-M(/11-/13) (/12-/14)(/13-/15) (/l4-M(/l5-/ll)' ^ ^ 

However, in this form, unlike the previous one, the terms attached to the variable k (say), contain 
not only the fields on two edges zq, ea parallel to the A;-axis. This form is therefore less suited for 
the solution of the Cauchy problem for the constrained f gh-system, which we discuss now. 



Theorem 17 For arbitrary a,/3 G C the constraint (5j_5), ( |5. (\ ) is compatible with the equations 



Proof. To prove this statement, one has to demonstrate the solvability of a reasonably posed 



Cauchy problem for the f gh-system constrained by (5.5), ( |5.6| ). In this context, it is unnatural to 
assume that the fields u, v vanish at the origin, so that we replace (only in this proof) the left-hand 
sides of ( |5.5D , ( |5.6| ) by au + cp, f3v + tp, with arbitrary 0, ■0 G C We show that reasonable Cauchy 
data are given by the values of two fields u, v (say) in three points 30) 3i = 3o + 1) 32 = 3o + ^, 
where 30 is arbitrary. According to Lemma & these data yield via the equations of the f gh-system 



the values of u, u in 33 = 30 + 1 + ^- Further, these data together with the constraint (|5.5|) , ( p.q ) 
determine uniquely the values of u, v in 34 = 30 + u;^. Indeed, assign n(34) = ^, f (34) = rj, where 
^, rj are two arbitrary complex numbers. The constraint uniquely defines the values of u, v in the 
point 35 =30 — 0;. The requirement that these values agree with the ones obtained via Lemma |8| 
from the points 30, 31, 34, gives us two equations for ^, r]. It is shown by a direct computation that 
these equations have a unique solution, which is expressed via rational functions of the data at 30, 
31, 32. It is also shown that the same solution is obtained, if we work with 35 = 3o ~ 1 instead of 
35. Having found the fields u, f at 34, we determine simultaneously u, f at 35, 35. Now a similar 
procedure allows us to determine u, t> at 37 = 30 + 2 and 33 = 30 + 2a;, using the constraint at the 
points 31 and 32, respectively. Simultaneously the values of u, v are found at 39 = 30 + 2 + cj and 
310 = 30 + 1 + 2a;. A continuation of this procedure delivers the values of u, v on the both semiaxes 



/c : /c > 0| U I3 = ^cj : ^ > 0|, 



using the condition that the constraint (|5.5| ), (|5.6| ) is fulfilled on these semiaxes. As we know from 
Proposition 0, these data are enough to determine the solution of the f gh-system in the whole 
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sector 



It remains to prove that this solution fulfills also the constraint (|57 
This follows by induction from the following statement: 



{3 = A; + ^w : > 0} = {3 € V{TC) : < arg(3) < 2^/3}. 

in the whole sector. 



Lemma 18 // the constraint (5^), ([5.(\ ) is satisfied in 30; di, I2, then it is satisfied also in 33. 



The constraint at 33 includes the data at five points 31, 32, 33, 39, 310- As we have seen, the data 
at 33 1 39 5 310 are certain (complicated) functions of the data at 30, 31, 32- Therefore, to check the 
constraint at 33, one has to check that two (complicated) equations for the values of u, v at 30, 31, 
32 are satisfied identically. This has been done with the help of the Mathematica computer algebra 
system. I 




Figure 4: To the proof of Theorem 17: labelling of the points 



Now we show how the constraint (|5.5| ), (|5.q ) appears in the context of isomonodromic solutions 
of integrable systems. In this context, the results look better with a different gauge of the transition 
matrices for the /g/i-system. Namely, we conjugate them with the matrix diag(l, A, A^), and then 
multiply by (1 + A^)^/^ in order to get rid of the normalization of the determinant. Writing then 
/i for A^, we end up with the matrices 




fgh = 1. 



The zero curvature condition turns into 

'C(e3,At)/3(e2,/i)>C(ei,/i) = (1 +fi)I, 



(5.13) 



(5.14) 



^1) ^2, ^3 being the consecutive positively oriented edges of an elementary triangle of TC. This 
implies some slight modifications also for the notion of the wave function. Namely, the previous 
formula does not allow to define the function ^ on V{TC) such that 

^(32,/") =>C(e,/x)^(3i,//) 

holds, whenever e = (31,32)- The way around this difficulty is the following. We define the wave 
function ^ on a covering of V(TC). Namely, over each point 3 = /c + + muj'^ now sits a sequence 



(/u) = (1 +^)"^fc,^,„(;u), nGZ. 



(5.15) 
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The values of these functions in neighboring vertices are related by natural formulas 



^ k+i,e,m{fJ') = >C(eo,/x)^'fc,^,m(/"), eo = {},d + 1), 

^fc,£+l,m(^) = >C(e2,/x)^'fc,^,m(/W), ^2 = {d,d + 
*fe/,m+l(^) = >C(e4,/i)^'A,.,f,m(/U), 64 = (3,3 + ^^). 



(5.16) 



We call a solution {u,v) : V{TC) i— > of the equations (^.lOP isomonodromic (cf. [^), if there 
exists the wave function ^ : Z-^ ^ GL(3, C)[/u] satisfying and some linear differential equation 

in 



(5.17) 



where Ak,e,m{fJ-) are 3x3 matrices, meromorphic in ^u, with the poles whose position and order do 
not depend on k, i, m. 



Obviously, due to ( 5.15 ), the matrix A has to fulfill the condition 

n 



Ak+n,e+n,m+n{^J') — Ak,^,m{^J') + 



n e z. 



(5.18) 



Theorem 19 Solutions of the equations ( 4-l(\ ) satisfying the constraints ^5.3i) , (5.6) are 
odromic. The corresponding matrix Ak tm ''■^ given by the following formulas: 



A, 



k,£,m 



Ck,e,m Djl) 



1 + /U /U 

where Ck/^m o-nd D{i) are ^-independent matrices: 
Pq 2 4 cire rank 1 matrices 



1 



fjdj + S'j/j+S + fj+39j+3 



fjdj -fj9jfj+3 fj9jfj+3gj+3 
-9j 9jfj+3 -9jf 3 +393 +3 
1 ~/j+3 f 3 +393 +3 



and the matrix D is well defined on V{TC) and not only on its covering 1? : 

m = 



-(2a + /3)/3 au f3a — aa' 

{a-(3)/3 (3v 

(2/3 + a)/3 



where the functions a, a' : V{TC) i-^ C are solutions of the equations ( \^.16 ), (4-li)- 
Proof can be found in the Appendix ^ I 



isomon- 



(5.19) 



(5.20) 



0,2,4, (5.21) 



(5.22) 



6 Isomonodromic solutions and circle patterns 

We now consider isomonodromic solutions of the fghsystem satisfying the constraint (|5.5| ), (^.fip, 
which are special in two respects: 

• First, the constants a and P in the constraint equations are not arbitrary, but are equal: 
a = (3, so that 7 = 1 — 2a. 
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Second, the initial conditions will be choosen in a special way. 



We will show that the resulting solutions lead to hexagonal circle patterns. 

First of all, we discuss the Cauchy data which allow one to determine a solution of the fgh- 



system augmented by the constraints (5.5), (p^). Of course, the fields u, v, w have to vanish in 
the origin 3 = 0. Next, one sees easily that, given u and v in one of the points neighboring to 
0, the constraint allows to calculate one after another the values of u and v in all points of the 
corresponding axis. For instance, fixing some values of n(l) and v{l), we can calculate all u{k) and 
v{k) from the relations 



au{k) = k 



f{k)g{k)f{k - 1) 



f{k)g{k) + g{k)f{k - 1) + f{k - l)g{k - 1) ' 



(6.1) 



I5v{k) = k 



where we have set 



g{k)f{k-l)g{k-l) 



f{k)g{k) + g{k)f{k - 1) + f{k - l)g{k - 1) ' 



f{k) = u{k + 1) - u{k), g{k) = v{k + 1) - t;(A;). 



(6.2) 



(6.3) 



Indeed, we start with u{0) = 0, v{0) = 0, /(O) = n(l), g{0) = v{l), and continue via the recurrent 
formulas, which are easily seen to be equivalent to (|6.l[) , (|6.2|) , (l6.3|) : 

u{k) =u{k-l) + f{k-l), v{k)=v{k-l)+g{k-l), (6.4) 



fik) 
9{k) 



au{k) 
(3v{k) 



I3v{k) 



au{k) (3v{k) 



(6.5) 
(6.6) 



k 



f{k-l) g{k-l) 



So, given the values of the fields u and v (and hence of w) in the points 3 = 1 and 3 = uj, we get 
their values in all points i = k and 3 = of the positive k- and ^-semiaxes. It is easy to see that 
u{k)/u{l) and v{k)/v{l) do not depend on n(l) and f(l), respectively, so that all points u{k) lie 
on a straight line, and so do all points v{k). Similar statements hold also for all points u{(.lo) and 
for all points v{£uj). And, of course, the third field w behaves analogously. 
So, we get the values of u and v in all points on the border of the sector 



= I3 G V{TC) : < arg(3) < 27r/3| = = k + lu; : k,l > Oy 



(6.7) 



Proposition ^ assures that these data determine the values of u and v in all points of S. By Theorem 
17 (more precisely, by Lemma 18) the solution thus obtained will satisfy the constraint (^), ( ^ ) 
in the whole sector S. 

Now we are in a position to specify the above mentioned isomonodromic solutions. 

Theorem 20 Let (3 = a. Let u,v,w : S ^ C he the solutions of the f gh-system with the constraint 



(5.5), (5.6), with the initial conditions 



ti(l) = v{l) = 1, u{uj) = v{uj) = exp{i9), 



(6.8) 



where < < it. Then all three maps u, v, w define hexagonal circle patterns with MR = —1 in the 
sector S. More precisely, if = 3' + s'^; k = 1,2, . . . ,6, are the vertices of an elementary hexagon 
in this sector, then: 
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• u{}i),u{^2), ■ ■ ■ ,u{}6) lie on a circle with the center in u{i') whenever 3' G 5 \ V{TLCi), 

• f (31), ^(32), . . . ,w(36) lie on a circle with the center in v[i') whenever i' € S\V{1-LC2), 

• ?y(3i), ^(32)1 • • • ,w{])^) lie on a circle with the center in w{i') whenever 3' € 5 \ V{'HCq). 

Proof follows from the above inductive construction with the help of two lemmas. The first 
one shows that \i [3 = a then the constraint yields a very special property of the sequences of the 
values of the fields u, v, w in the points of the k- and £-axes. 

Lemma 21 If (3 = a, then for k,i > I: 

\u{3k - 1) - u{3k - 2)\ = \u{3k - 2) - u{3k - 3)\, (6.9) 
\vi3k) - v{3k - 1)\ = \vi3k - 1) - v{3k - 2)\, (6.10) 
\w{3k + I) - w{3k)\ = |'u;(3A;) -'^(SA; - 1)1, (6.11) 

\u{{3i - l)uj) - u{{3i - 2)uj)\ = \u{{3i - 2)uj) - u{{3i - 3)u;)\, (6.12) 
\v{3eLo) - v{{3i - l)uj)\ = \v{{3e - l)uj) - v{{3i - 2)lo)\, (6.13) 
\w{{3i + l)io) - w{3iuj)\ = \w{3iuj) - w{{3i - l)uj)\. (6.14) 

The second one allows to extend inductively these special properties to the whole sector (|6.7|). 

Lemma 22 Consider two elementary triangles with the vertices 30, 31 = 3o + 1; 32 = 3o + o-nd 
33 = 30 + 1 + i^- Suppose that 

(i) \u{li) - u(3o)| = 1^(32) - ^i(3o)|; 

(ii) (31)^(30)^^(32) = and l.u{])i)u{])o)u{i2) = 2tt -2'd for some ■!?. 
Then 

^(33) - u{iq)\ = |u(3i) - u(3o)| = \u{i2) - u{io)\, (6.15) 

and hence 

^(33) - ^^(3i)| = 1^^(30) - i^(3i)|, ^(33) - -^(32)1 = |w(3o) - ^(32)1 (6.16) 

and 

^(33) - ^i'(3i)l = ^(33) - ^^(32)1 = ^(33) - w(3o)| (6.17) 

The assertion of this lemma is illustrated on Fig. ^. 

First of all, we show how do these lemmas work towards the proof of Theorem 20. The initial 
conditions (|6.8|) imply: 



w{l) = 1, w{uj) = exp(-2i6') = exp(i(27r - 261)). (6.18) 

Therefore, the conditions of Lemma |2^ are fulfilled in the point 30 = with the fields (u), u, v) 
instead of {u,v,w). From this Lemma it follows that 

(ao) The points w{l), w{u)), w{l + uj) are equidistant from w{0); 

(bo) The points v{0), v{l), v{uj) are equidistant from v{l + uj); 
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(cq) The points u(l +uj), u{0) are equidistant from u{uj); 
(do) The points u(l +uj), u{0) are equidistant from n(l). 

Since, by Lemma ^ we have \u{0) — u{l)\ = \u{2) — u{l)\, there fohows from (do) that |u(l +uj) — 
u{l)\ = \u{2) — 14(1)1. Finally, from Lemma there follows that (see Fig. |^ 

Zv{2)v{l)v{l + w) = TT - Vi, lu{2)u{l)u{l + oj) = TT - (pi = 2'ipi = 2tt - 2{tt - V'l). 

Therefore the conditions of Lemma are fulfilled in the point 30 = 1 with the fields {u,v,w). We 
deduce that 

(ai) The points u{2), u{l + u>), u{2 + w) are equidistant from u{l); 
(bi) The points w{l), w{2), w{l + lo) are equidistant from w{2 + to); 

(ci) The points v{2 + lu), v{1) are equidistant from v{l + w), which adds the point v{2 + uj) to the 
list of equidistant neighbors of v{l + uj) from the conclusion (bo) above; and 

(di) The points v{2 + w), v{l) are equidistant from v{2). 

By Lemma |2l|, we have \v{l) — v{2)\ = \v{3) — v{2)\, and there follows from (di) that \v{2 + uj) — 
v{2)\ = \v{3) — v{2)\. Finally, from Lemma ^ there follows that (see Fig. ^) 

Zw{3)w{2)w{2 + w) = TT - Vs, lv{3)v{2)v{2 + w) = vr - </>3 = 2^3 = 27r - 2(7r - Vs). 



Hence, the conditions of Lemma 22 are again fulfilled in the point 30 = 2 with the fields {v,w,u). 

These arguments may be continued by induction along the fc-axis, and, by symmetry, along the 
-^-axis. This delivers all the necessary relations which involve the points ^ = k + iuj with /c < 1 or 
i < 1. We call them the relations of the level 1. 

The arguments of the level 2 start with the pair of fields {v, w) at the point 3 = 1 + 0;. We have 
the level 1 relation 

\v{2 + uj) - v{l + uj)\ = \v{l + 2uj) -v{l+uj)\. 

For the angles, we have from the level 1 (see Fig. |6|): 

Zw{2 + uj)w{l + uj)w{l + 2uj) = 27r - (^1 + + V'4 + V's), 
lv{2+uj)v{l+uj)v{l+2uj) = 27r - (01 +02 + </'4 + </'5) = 27r - 2(27r- ^1 - 7/^2 - ^^4 - V's)- 
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Figure 6: To the proof of Theorem 2C: similar isosceles triangles for u, v, and w 
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So, the conditions of Lemma 22 are again satisfied in the point = 1 + u for the fields {v,w,u). 
Continuing this sort of arguments, we prove ah the necessary relations which involve the points 
^ = k + iio with k < 2 or i < 2, and which will be called the relations of the level 2. The induction 
with respect to the level finishes the proof. I 



It remains to prove Lemmas 21 and 22 above 



What concerns the key Lemma |2^, it might be instructive to give two proofs for it, an analytic 
and a geometric ones. The first one is shorter, but the second one seems to provide more insight 
into the geometry. 



Analytic proof of Lemma 22. We rewrite the assumptions of the lemma as 



U2 - Uq 



(ni-no)e2^(--'') = (ni-txo)e- 



■211) 



and 



V2-vo = c{vi - vo)e 



c> 0. 



Geometric proof of Lemma 22. The equations of the /g'/i-system imply that the triangles 
uqUiu^ and viV2,vq are similar, and the triangles uqU2U^ and V2V'iVQ are similar. Therefore, 



\Vi - Vo\ 



\Vi - fsl 



From |uo — nil 



\uq-U2,\ \uo-ui\ 
U2\ there follows now 

1^1 - Wo I 



\Vl - Vi\ 



\V2 - vo\ 



\V2 - V0\ 
\V2 - 



\V2 - VZ\ 
\uq - U2\ 



(6.19) 




Figure 7: To the proof of Lemma 



^^1 



Denoting the angles as on Fig. 0, we have: 

XI + X2 = ^9, 01 + <^2 = 27r - 2^9, 

hence 

^1 + V'2 = 27r - ((/>i + 02) - (Xl +X2)='& = Xl+ X2- 

In other words, 

ZV1V3V2 = ZV1V0V2. 



(6.20) 
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The relations ( |6.19| ), ( |6.20| ) yield that the triangles viVsV2 and viVoV2 are similar. But they have 
a common edge [^1,^2], therefore they are congruent (symmetric with respect to this edge). This 
implies that the triangles V0V2V3 and vqViVs are isosceles, so that xi = "01 and X2 = and 

\vo -Vi\ = \V3 -Vi\, \vo -V2\ = \vs - V2\. 

Therefore 

1^3 - '"ol = \U1 - Uo\ = \U2 - Uo\- 

Lemma is proved. I 

As for Lemma its statement is a small part of the following theorem and its corollary. 



Theorem 23 If (3 = a, then the recurrent relations i\6.^ , {6.5), (6.C) with n(l) = v{l) = 1 can 
be solved for u{k) , v{k), f{k), g{k) {k > 0) in a closed form: 

2k 2k + 2a 

= ——Ui{k), «(3fc + l) = — — -ni(fc), tt(3A: + 2) = 2ni(A:), (6.21) 
k + 2a k + 2a 

fiSk - 1) = f{3k) = f{3k + 1) = Uiik), (6.22) 

k + 2a 



and 



v{3k-l) = ^^U2{k), v{3k) = -^U2{k), v{3k + l) = U2{k), (6.23) 
k + a k + a 

g{3k - 2) = g{3k - 1) = g{3k) = U2ik), (6.24) 

k + a 



where 



^ (l + 2o)(2 + 2„).. (;c + 2a) ^ (1 + „)(2 + o) (^c + g) 

^ ^ {l-a){2-a)...{k-a) ^ ' {I - 2a){2 - 2a) . . . {k - 2a) ^ ' 

Proof. Elementary calculations show that the expressions above satisfy the recurrent relations 



(6.4), (p^), (|6.6D with (5 = a, as well as the initial conditions. The uniqueness of the solution 



yields the statement. We remark that similar formulas can be found also in the general case a ^ /?, 



however, the property formulated in Lemma 21 fails to hold in general 



Corollary 24 If (3 = a, and n(l) = v{l) = \, then for the third field w{k), h{k) {k > 0) we have: 

k — l + 2a„,,^ k ^ A; + l — 2q;„,,, 

^^ 3fc-l = ^ Hs , w{3k) = ——Il3{k), w{3k + l) = ^- Hsfe, 

1 — 2a 1 — 2a 1 — 2a 

(6.26) 

t- _l_ 1 —2a 

h{3k-l) = h{3k)=Il3{k), h{3k + l) = —- n3(A;), (6.27) 

k + a 

where 

^ (l-a)(2-a)...(fc-a) (1 - 2a)(2 - 2a) . . . (A: - 2a) 
' a{l + a)...{k-l + a) 2a(l + 2a) . . . (fc - 1 + 2a) ' ^ ' 

Proof. The formulas for h{k) = {f{k)g{k))-^ follow from (|]2|), d^H). The formulas for w{k) = 
w{k — 1) + h{k — 1) with w{Q) = follow by induction. I 
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7 Discrete hexagonal and log 2; 

Although the construction of the previous section always delivers hexagonal circle patterns with 
MR = —1, these do not always behave regularly. As a rule, they are not embedded (i.e. some 
elementary triangles overlap), and even not immersed (i.e. some neighboring triangles overlap), cf. 



Fig. |8|). However, there exists a choice of the initial values (i.e. of 9 in Theorem 2C) which assures 
that this is not the case. 




Figure 8: A non-immersed pattern with 6 ^ lira. 



Definition 25 Let < a = l3 < ^, so that 0<7 = 1 — 2a<l. Set 9 = 27ra. Then the hexagonal 
circle patterns of Theorem are called: 

u, V : the hexagonal z'^" with an intersection point at the origin; 
w : the hexagonal z^^ with a circle at the origin. 

In other words, for the hexagonal z'^" the opening angle of the image of the sector (6.7) is equal to 
27ra, exactly as for the analytic function z i— s- z^". 

Conjecture 26 For < a < ^ the hexagonal circle patterns z^'^ with an intersection point at the 
origin and z^^ with a circle at the origin are embedded. 

For the proof of a similar statement for circle patterns with the combinatorics of the square grid 
see [AE], where it is proven that they are immersed. 

Remark. Actually, the u and v versions of the hexagonal z^" with an intersection point at 
the origin are not essentially different. Indeed, it is not difficult to see that the half-sector of the 
u pattern, corresponding to < arg(3) < vr/S, being rotated by vra, coincides with the half-sector 
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of the V pattern, corresponding to vr/S < arg(3) < 27r/3, and vice versa. For the w pattern, both 
sectors are identical (up to the rotation by ttj). So, for every < a < ^ we have two essentially 
different hexagonal pattrens z^". 

It is important to notice the peculiarity of the case when a = n/N with n, S N. Then one can 
attach to the u, u-images of the sector S its copies, rotated each time by the angle 27ra = 2Trn/N. 
The resulting object will satisfy the conditions for the hexagonal circle pattern everywhere except 
the origin 3 = 0, which will be an intersection point of M = nN circles. Similarly, if 7/2 = n' /N', 
and we attach to the if-image of the sector S its N' copies, rotated each time by the angle 
27r7 = Ann' /N' , then the origin 3 = will be the center of a circle intersecting with M' = n'N' 
neighboring circles. See Fig. ^ for the examples of the tD-pattern with 7 = 1/5 and the n-pattern 
with Q = 1/5. 

Now we turn our attention to the limiting cases a = 1/2 and a = 0. 
7.1 Case a = |, 7 = 0: hexagonal z^^'^ and logz 

It is easy to see that the quantities g{k), k > 1, and v{k), k > 2, become singular as a — > ^ (see 
(6.24) and ( |6.23 )). As a compensation, the quantities h{k), k > 1, vanish with a — > ^, so that 



w{k) — > w{l) = 1 for all k > 2. Similar effects hold for the ^-axis, where v{iuj), i > 2, become 
singular, and w{iuj) 1 for all £ > 1. (Recall that for the w pattern we have: w{u>) = e^'^*'^ — > 1). 
These observations suggest the following rescaling: 

u = u, V = v/{l -2a), w = l + {l-2a)w. (7.1) 

In order to be able to go to the limit q — > ^, we have to calculate the values of our fields in several 
lattice points next to 3 = 0. Applying formulas ( [4.13 ), ( 4.11| ), we find: 



u(0) = 0, u{l) = l, u(u;) = e2^*°, u(l + w) = 1 + e^™, (7.2) 

g27riQ: 

1 + e 

W{0) =0, W{1) = 1, w{lu)= e2-^(l-2a) ^ + ^) ^ e-i{l-2a) ^7 4) 



^(0)=0, vil) = l, T;(a;) = e2™, ^(l+^) = __^, (7.3) 



For the rescaled variables n, v, w in the limit q — > ^ we find: 



u{0) = 0, u{l) = 1, u{uj) = -1, u{l + w) = 0, (7.5) 



v{0) = 0, v{l) = 0, v{uj) = 0, v{l + uj) = -, (7.6) 



vr 



w{0) = oo, 'w{l)=0, w{ijj) = 2'iTi, w{l + ijj) = Tii. (7-7) 

These initial values have to be supplemented by the values in all further points of the k- and £-axes. 
From the formulas of Theorem |2^ there follows: 



')ki,\ 9^k^ 9^k^ 

^^'^^ = (2^^ • ^^'^ + '^=(2fc3I)n-(2^ + l)' ^(3A: + 2) = ^^^-^.(2A: + 2), 



(7.8) 



/(3fc - 1) = f{3k) = }{?,k + 1) = j^^, (7.9) 
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Figure 9: The hexagonal patterns 7?!^ with a circle at the origin and with an intersection point at 
the origin. 
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Figure 10: Some examples of u;-pattern: 7 = 1, 2/3, 1/2, 2/5, 1/3, 2/7. 
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and 

(7.10) 

giSk - 2) = g{3k - 1) = g{3k) = (7.11) 



which have to be augmented by u{ku;) = —u{k), v{kuj) = —v{k). From Corohary ^ there follow 
the formulas for the edges of the w lattice: 

h(2,k - 1) = h{3k) = h({3k - l)u;) = hi3kuj) = \, k>l, (7.12) 

k 

h{3k + l) = h{{3k + l)uj) = \ fc>0. (7.13) 

Definition 27 The hexagonal circle patterns corresponding to the solutions of the fgh-system in 
the sector (\6. defined by the boundary values ( 7.5)-( '7A^ ) are called: 



u, V : the hexagonal z^l"^ with an intersection point at the origin; 
w : the symmetric hexagonal logz. 



o o o 



Alternatively, one could define the lattices u, v, w as the solutions of the f gh-system with the 
initial valu es ([7.5D-(|7.7] ), satisfying the constraint (|5.5D, ( |5.6| ) with a = P = 1/2. In this appoach 
the values (|7.8D -( [r.l3D would be derived from the constraint. Notice also that the formulas (^3), 
(|5.9| ) in this case turns into 

1 = k— — -i — ^ + — — ^ + — — ^ (^-^^^ 

fo9o + 9of3 + fsds /252 + 52/5 + fbh hi A + 54/1 + 

000 000 000 

o o 00 00 o o 00 00 "'""'■o o 00 o o • \'-^^J 

^0/0 + /o^3 + /13/3 /12/2 + /2/15 + ^^5/5 ^4/4 + /4/1I + ^1/1 

7.2 Case a = 0, 7 = 1: hexagonal log 2; and 

Considerations similar to those of the previous subsection show that, as a — > 0, the quantities h{k), 
k > 1, and w{k), k > 2, become singular (see (6.27) and ( 6.26| )). As a compensation, the quantities 



f{k), k > 2, and g{k), k > 1, vanish with a ^ 0, so that u{k) u{2) = 2 for all A; > 3, and 
v{k) — > v{l) = 1 for all k > 2. Similar effects hold for the ^-axis. These observations suggest the 
following rescaling: 

u = 2 + 2au, v = l + av, w = w/{2a^). (7.16) 

It turns out that in this case we need to calculate the values of these functions in a larger number 
of lattice points in the vicinity of 3 = 0. To this end, we add to ([7^)-([7r^) the following values, 
which are obtained by a direct calculation: 

1 -I- p27rio 

u{2) = 2, n(2^) = 2e2™, u{2 + u:) = ^ ^ J^^,^,^ _ ,y (7.17) 
+ + "P + ^)=i32S(l + e-"), (7.18) 
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Figure 11: The patterns z^/^ with an intersection point at the origin, and the symmetric hexagonal 
logz; the second pattern coincides with the first one upon the rotation by 7r/2 
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vi2) 



1 — a 



2a' 
v{l + 2uj) 



v{2uj) 



1 



1 -2a 

„27ria 



1 + a(e2™ - 1)' 



"^' + "^^= l + a(e-'--l) ' 

2g27r'ia 

t;(2 + 2lo) = 



1 + e^'T*" 



w;(2) 



1 — a 



a 



■u;(2w) 



w{l + 2w) 



1 — a 



a 



-27rja 



a(e- 



1)^ 



w{2 + uj) 



w{2 + 2a;) 



27r«a 



a(e 
1 — a 



-27ria 



a 



(7.19) 
(7.20) 

(7.21) 
(7.22) 



From ( [7.2D ~( [7r^ ) and ( 7.17| )-( [7.22| ) we obtain in the limit a — > under the rescahng (7.16) the 
following initial values: 

o,> o,> o,, o,> o,, 

n(0) = oo, n(l) = oo, uiuj) = oo, u{2) = 0, u[2ui) = 2m, (7.23) 
u{l+uj) = m, u{2 + uj) =7ri, u{l + 2uj) = ni, u{2 + 2uj) = 1 + ni, (7.24) 



v{0)=oo, v{l) = 0, v{uj) = 2-Ki, v{2) = l, t;(2w) = 1 + 27ri, 
v{l + uj) = oo, v{2 + uj) = {), v{l + 2Lo) = 2'Ki, v{2 + 2Lo)=m, 

^(0) = 0, w{l) = 0, w{uj) = 0, w{2) = 0, w{2uj) = 0, 



w(l + w) = 0, 'u;(2 + w) = -, ti;(l + 2u;) 



TT 



t(;(2 + 2a;) = 0. 



(7.25) 
(7.26) 

(7.27) 
(7.28) 



These initial values have to be supplemented by the values in all further points of the k- and ^-axes. 
From the formulas of Theorem |2^ there follow the expressions for the edges of the lattices u, v: 



fi3k-l) = fi3k) = f{3k + l) = f{{3k-l)u) = f{3kuo) = f{{3k + l)Lo) 
g{3k - 2) = °g{3k - 1) = g{3k) = °g{{3k - 2)uj) = g{{3k - l)a;) = ^(S/ca;) 



1 

k' 
1 



The formulas of Corollary |24| yield the results for the lattice w: 

w{3k) = k^, w{3k + l) = k'^{k + l), w{3k + 2) = k{k + if , k>l, 

so that 

^(3A; - 1) = ^(3A;) = A;2, h{3k + 1) = k{k + 1), k>l. 
Of course, one has also w{kuj) = w{k). 

Definition 28 The hexagonal circle patterns corresponding to the solutions of the fgh-system in 
the sector ( \6. ?[ ) defined by the boundary values ( 7.2c )-( 7.33( ) are called: 



k>l, 

(7.29) 
k>l. 

(7.30) 

(7.31) 
(7.32) 



u, V : the asymmetric hexagonal log z; 
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w : the hexagonal with a (degenerate) circle at the origin. 

It is meant that the -u-image of the half-sector < arg(3) < tt/S is not symmetric with respect 
to the hne ^(li) = 7ri/2 (the image of arg(3) = vr/6, and the same for v. Instead, this symmetry 
interchanges the u pattern and the v pattern, see Fig. [l^ . 

Alternatively, one can define these lattices as the solutions of the fgh-system with the initial 
values ( [7.23|) - ([?'.28 ), satisfying the constraint (^^), (|5.6D , which in the present situation degenerates 
into 

°o° °o° °o° 

~ '°~i> 

fo9o + 50/3 + /353 f292 + 52/5 + /sffs h94 + 54/1 + fl9l 
1 = „ 'ff' o ^ ^ 'i'h o ^ ^ ''/l'' „ ^ ■ (7.34) 

/o5o + 50/3 + f393 f292 + 52/5 + f5°95 hi A + Ijl + fill 

Just as in the non-degenerate case, these formulas allow one to calculate inductively the values of 
tt, V on the k- and £-axes. The formulas ( ^.71) , (^.91) hold literally with 7=1. 



8 Conclusions 

In this paper we introduced the notion of hexagonal circle patterns, and studied in some detail a 
subclass consisting of circle patterns with the property that six intersection points on each circle 
have the multi-ratio —1. We established the connection of this subclass with integrable systems 
on the regular triangular lattice, and used this connection to describe some Backlund-like trans- 



formations of hexagonal circle patterns (transformation u ^ v ^ w, see Theorems 13, 14), and to 
find discrete analogs of the functions z", logz. Of course, this is only the beginning of the story 
of hexagonal circle patterns. In a subsequent publication we shall demonstrate that there exists 
another subclass related to integrable systems, namely the patterns with fixed intersection angles. 
The intersection of both subclasses constitute conformally symmetric patterns, including analogs 
of Doyle's spirals (cf. pH| ). 

A very interesting question is, what part of the theory of integrable circle patterns can be 
applied to hexagonal circle packings. This also will be a subject of our investigation. 

This research was financially supported by DFG (Sonderforschungsbereich 288 "Differential 
Geometry and Quantum Physics"). 



A Appendix: Square lattice version of the /(//i— system 

Dropping all edges of E{TC) parallel to the m-axis, we end up with the cell complex isomorphic 
to the regular square lattice: its vertices i = k + lui may be identified with {k,€} G I?, its edges 
are then identified with those pairs [(A;i,^i), (A;2,^2)] for which l^i — /c2[ + |^i — ^2! = 1, and its 
2-cells (parallelograms) are identified with the elementary squares of the square lattice. Hence, flat 
connections on TC form a subclass of fiat connections on the square lattice. A natural question is, 
whether this inclusion is strict, i.e. whether there exist flat connections on the square lattice which 
cannot be extended to fiat connections on TC. At least for the fgh-system, the answer is negative: 
denote hy M- C SL(3,C)[A] the set of matrices ( |4.2| ), then fiat connections on the regular square 
grid with values in M are essentially in a one-to-one correspondence with fiat connections on TC 
with values in Ai, i.e. with solutions of the fgh-system. This is a consequence of the following 
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Figure 12: The asymmetric patterns log z and the hexagonal pattern ^ with a circle at the origin; 
the upper half of the first pattern coincides with the lower half of the second one, and vice versa 
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statement dealing with an elementary square of the regular square lattice: a flat connection on 
such an elementary square with values in M can be extended by an element of A4 sitting on its 
diagonal without violating the flatness property. More precisely: 

Lemma 29 Let 



L1L2 = I/3L4, where Li e M (i = 1,2,3,4), 



and let the off-diagonal parts of Li, L2 be componentwise distinct from the off-diagonal parts of 
L3, L4, respectively. Then there exists Lq ^ M such that 



L0L1L2 — LqLsL^ 



L. 



La 




Lo 



Figure 13: To Lemma ^ 



Proof. We have to prove that (L1L2) ^ = {L3L4) ^ G A4. It is easy to see that it is necessary 
and sufficient to prove that the entries 13, 21, 32 of this matrix vanish, i.e. that there holds 



fi9i + /2C/1 + /252 = /353 + fm + fm = 0, 



(A.l) 



as well as two similar equations resulting by two successive permutations {f,g,h) 1— > {g,h,f). We 
are given the relations figihi = 1 and 



/i + /2 = /a + /4, 91+92= 91,+ 94, hi + h2 = h3 + 



(A.2) 



fl92 = /354, 9lh2 = 53^4, /il/2 = hsf4. 

In order to prove (A.l), we start with the third equation in (|A.2|): 



/ii 1 



hs 
hi 



Using figihi = 1 and ( |A.3|) , we find: 

hs ^ /igi ^ g4gi fu ^ gi_ 
hi fz9z 929?,' h2 93 



Plugging this into ( A. 4 ), we get: 



g2g3 - glg4 _ 91 - 93 
fl9l9293 f292g3 ' 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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Now, due to the second equation in ( [A. 21) , we find: 



5253 - 9194 = 92(93 - 91) + 91(92 - 9a) = (91 + 92)(g3 - 9i)- 
Substituting this into (A. 6), we come to the equation: 



/ ^ I 91+92 I 



0. 



(A.7) 



(A. 



Since, by condition, gi ^ (73, we obtain f2(gi + 92) + fi9i = 0, which is the equation ( A.l) 



This result shows that the /g/i-system could be alternatively studied in a more common frame- 
work of integrable systems on a square lattice. However, such an approach would hide a rich and 
interesting geometric structures immanently connected with the triangular lattice. It should be said 
at this point that the one-field equation ( [4.19 ) was first found, under the name of the "Schwarzian 
lattice Bussinesq equation" by Nijhoff in Q using a (different) Lax representation on the square 
lattice. The same holds for the one-field form of the constraint ( |5.10 ). 



B Appendix: Proofs of statements of Sect. ^ 



Proof of Proposition 15. The arguments are similar for both equations (5^), (|5.6| ). For instance, 
for the first one we have to demonstrate that 



fo9oh 



+ 



7292/5 



fo9o + 50/3 + f393 ' f292 + 52/5 + f595 ^ f49A + 9ifl + fl9l 
/0/3 _| 72/5 |_ /4/1 



fmfi 



0. 



(B.l) 



/O + /3 + /353/5O f2 + h + h9b/92 fA + h+ fl9l/9A 

To eliminate the fields g from this equation, consider six elementary triangles surrounding the 
vertex 3. The equations (4^) imply: 



91 _ 


fo + h 


92 _ 


fi 


93 _ 


/2 + /3 


90 


fi ' 


91 


/1+/2' 


92 


/3 


95 _ 


/5 + /0 


94 _ 


/5 


93 _ 


/3+/4 


90 


h ' 


95 


/4 + /5' 


92 


/3 



Therefore, 



/O + /3 + h 



93 

s — 
90 



/0 + /3 
/0 + /3 



(/0+/l)(/2 + /3) (/0-/2)(/l-/3) 



/1 + /2 



/1 + /2 



(/5 + /o)(/3 + /4) (/4-/o)(/3-/5 



/4 + /5 



(B.2) 
(B.3) 



/4 + /5 

By the way, this again yields the property MR = —1 of the lattice u, which can be written now as 

(fo + /l)(/2 + /3)(/4 + /s) = (fl + /2)(/3 + U)(f5 + /o), (B.4) 



Using ( [B.2D , an analogous expression along the £~axis, and an expression analogous to (|B.3| ) along 
the m-axis, we rewrite (IB. 11) as 



/o/3(/l + /2) 



+ 



/2/5(/3 + /4) 



+ 



/4/l(/2 + /3 



(/0-/2)(/l-/3) (/2-/4)(/3-/5) (/2-/4)(/l-/3) 



0. 



(B.5) 
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Clearing denominators, we put it in the equivalent form 

/o/3(/l + /2)(/2 - /4)(/3 - h) + /2/5(/3 + U){fo " /2)(/l " /s) 

+/4/l(/2 + /3)(/0-/2)(/3-/5) = 0. 

But the polynomial on the left-hand side of the last formula is equal to 

/2/3((/l + /2)(/3 + U){h + fo) - (/O + /l)(/2 + /3)(/4 + h 



and hence vanishes in virtue of ( [B.4D . I 

Proof of Proposition Denote the right-hand sides of (^), ( |5.6D , (|5.7D through U{^), 
V{^), W{i), respectively. In order to prove ( |5.7| ), i.e. = W{i), it is necessary and sufficient to 
demonstrate that 

7/10 = W(j) - Wii), 7/12 = W{i) - W{i), 7/14 = W{i) - W{i), 

(or, actually, any two of these three equations). We perform the proof for the first one only, since 
for the other two everything is similar. In dealing with our constraints we are free to choose any 
representative {k,i,m) for 3. In order to keep things shorter, we always assume in this proof that 
m = 0. Writing the formula 

in long hand, we have to prove that 

j = l-a-P = (k + l) ^^ ^[^2 -k- 



fm + gofs + hgz fo9o + 50/3 + fsgs 

'-^i^^^-i , '^)\, . (B.6) 

/252 + 52/5 + /555 J ^92 + 5-2/5 + /sffS 



Assuming that (5.5) and (|5.6|) hold, we have: 



a 



/o V ^ go 

Taking into account that fs = fo, gs = go, we find: 



.a n , T\ jyjuyi + goh , 50/3 + 7353 

a + p = (A; + i)^^ — k 



+1 

or, equivalently. 



fogo + goh + hgs fo9o + goh + fsgs 

( 7252 h/fo+ g2 /sffS / go _ /2g2/5//o +g2/5g5/g0 

V 7252 + 52/5 + 7555 hg2 + g2h + hg5 



1 o /; , i\ /353 , fogo 

7 = 1 — a — p = (fe + l)- 



/o5o + 50/3 + /353 /offo + 50/3 + /353 
^ / 7252 hlfo+ 52 755 / go _ /252/5/ /o + 52/555/50 \ ,^ 
V 752 + 527 + 755 h92 + 52/5 + /555 I ■ 



36 



The first two terms on the right-hand side already have the required form, since f^g^ = fogo = 1/hQ. 
So, it remains to prove that 



1292/5/ fo + 92/595/90 _^ 1292/5/ fo + 92/595/90 



1//10 



/292 +92/5 + /595 



/292 + 92/5 + /595 



/292 + 92/5 + /595 /292 + 92/5 + /595 

(Bi 



The most direct and unambiguous way to do this is to notice that everything here may be expressed 
with the help of the /(7/1-equations in terms of a single field h. After straightforward calculations 
one obtains: 



7~ /5 , 7~ 92 
/292-r + /595— = 
JO 90 


1 

/15 


1 ho{ho - /15) 


(B.9) 


f /s , J. 92 
/292-r + /595— = 

Jo 90 


1 


ho{ho - h2) 

\ — ; 


(B.IO) 


/292 +92/5 + /555 = 


{ho- 


- h5){hi - h2) 


(B.ll) 






/292 +92/5 + /595 = 


{ho- 


- h2){hi - /is) 


(B.12) 




/12/12/15 



Taking into account that /14 — /i2 = /ii — /is, we see that (B^) and Proposition 16 are proved. 



Proof of Theorem |l^. In order for the isomonodromy property to hold, the following com- 
patibility conditions of ( 5.1(]| ) with (5.17) are necessary and sufficient: ( ^.14| ) and 



^(^0) = A+i/,m>C(eo, /u) - C{to,ix)Ak. 



d_ 

— £(e2,^) = Ak,l+l,mC,{'^2-, - C{Z2, H)Ak,l,m, 



(B.13) 



Substituting the ansatz ( ^.19|) and calculating the residues at /i = —1, = and 11 = 00, we see 
that the above system is equivalent to the following nine matrix equations: 



Ck+l/^m^i'^O', — 1) 
Cfc,^+l,m'^(f2; — 1) 
Ck/,m+lC{u, — 1) 



C{zA,-l)Ck,t 



(B.14) 
(B.15) 
(B.16) 



^(3}^(eo,0) 
£'(T)^(e2,0) 
D{i)C{ti,{)) 



>C(eo,0)Z)(3), 
>C(e2,0)Z)(3), 
^(e4,0)D(3), 



(B.17) 
(B.18) 
(B.19) 



Ck+l,i,m + D(J))Q-Q{Ck,t,m + D{l)) = Q, 
Ck/+l,n^ + D(S)^Q-Q(Ck/,ra + D{l)^ = Q, 
Ck,i,m+l+D{^))Q-Q(Ck,i,,n + D{i)) = Q, 



(B.20) 
(B.21) 
(B.22) 
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where 



Q 




(B.23) 



We do not aim at solving these equations completely, but rather at finding a certain solution leading 
to the constraint (^), ( ^.61 ). The subsequent reasoning will be divided into several steps. 

Step 1. Consistency of the ansatz for Ck/^m- First of all, we have to convince ourselves 
that the ansatz ( ^.20 ), ( 5.21| ) does not violate the necessary condition ( 5.1§| ), i.e. that 

P2 + Pi + PQ = I. (B.24) 

Notice that the entries 12 and 23 of this matrix equation are nothing but the content of Proposition 



15. Upon the cyclic permutation of the fields {f,g,h) i— > {g,h,f) this gives also the entry 31. To 



check the entry 21, we proceed as in the proof of Proposition |15|. We have to prove that 
50 , 92 , 94 



+ 



+ 



fodO + 50/3 + /353 /252 + 52/5 + /555 fA9A + difl + fl9l 

1 1 1 

/o + /s + /353/90 f2 + h + h9^/92 fA + h + fl9l/9A 

/1 + /2 , h + h . f2 + h 



+ 



+ 



0. 



(/O - /2)(/l - /3) (/2 - h){h - h) {h - k){h - h) 

Clearing denominators, we put it in the equivalent form 

(/l + /2)(/2 - k){h - h) + (/3 + /4)(/0 - /2)(/l - /3) + (/2 + /3)(/o " /2)(/3 " h) = 0- 

But the polynomial on the left-hand side is equal to 

(/i + /2)(/3 + h){h + h) - (/o + /i)(/2 + h){h + h), 



and vanishes due to (B^). Via the cyclic permutation of fields this proves also the entries 32 and 
13 of the matrix identity ( 6.24] ). Finally, turning to the diagonal entries, we consider, for the sake 
of definiteness, the entry 22. We have to prove that 



/350 



+ 



f592 



+ 



fl9A 



fo9o + 9oh + f393 f292 + 52/5 + h95 f494 + 54/1 + fl9l 
h , f5 , fi 



+ 



+ 



/o + /s + /393/90 h + h + h9^/92 f4 + fi + fi9i/94 



/3(/l+/2) 



+ 



Mf3 + h) 



+ 



/l(/2 + /3) 



or 



(/0-/2)(/l-/3) {h-h){h-h) {h-h){h-h) 

Wl + /2)(/2 - /4)(/3 - h) + Uh + h){h - /2)(/l - /3) + 
/l(/2 + /3)(/0 - /2)(/3 - h) - (/O - /2)(/l - /3)(/2 " /4)(/3 " k) = 

Again, the polynomial on the left-hand side is equal to 

/3((/l + /2)(/3 + /4)(/5 + /o) - (/O + /l)(/2 + /3)(/4 + h)) , 
and vanishes due to ( [B.4| ). The formula (B.24) is proved. 
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Step 2. Checking the equations for the matrix Ck/^m- Next, we have to show that the 
ansatz (|5.2q ), ( ^.21| ) verifies ( p.l4| )- pl6| ). Notice that the matrices 




are degenerate, and that 



fg 



1 



and r] 



-f, fg 



are the right null-vector and the left null-vector of C{e, —1), respectively. In terms of these vectors 
one can write the projectors -Po,2,4 as 



Pi 



1 



j = 0,2,4.. 



Therefore we have: 



Poir)C{to,-l)=C{to,-l)Poii) = 0, 

P2(T)/:(e2,-i) =/:(e2, -1)^2(3) = 0, 

P4(3)>C(e4,-l) =/:(e4, -1)^4(3) = 0. 
In order to demonstrate (|B.14D -( [B.16| ) it is sufficient to prove that 

P2(J)C{zo, -1) = C{co, -1)^2(3), P4(^}C{co, -1) = C{co, 
Pi{7)C{c2, -1) = C{c2, -1)^4(3), i^o(3)^(e2, -1) = C{C2, 

Po(3)>C(e4, -1) = C{C4, -l)Po(3), i^2(3)>C(e4, -1) = C{c4, 



-1)A(3) = 0, 
-l)^'o(3) = 0, 

-1)^2(3) = 0. 



(B.25) 
(B.26) 
(B.27) 

(B.28) 
(B.29) 
(B.30) 



All these equations are verified in a similar manner, therefore we restrict ourselves to the first one. 



1 



1 



'C(eo, -'i-)C2vI 



or, in long hand, 
1 

/252 + 52/5 + hg5 



/2£2 
-52 
1 



^ I - ho/h^ 
_/o - fb 

\ h{gb - go) 



7252+52/5 + /595 



(/2 - fo)g2 
go - 52 

1 - /lo//l2 




To prove this we have, first, to check that these two rank one matrices are proportional, and then 
to check that their entries 31 (say) coincide. The second of these claims reads: 



1 - /lo//l5 



1 - ho/h2 



/252 + 52/5 + /555 /252 + 52/5 + /555 ' 

and follows from ( B.ll ), ( B.12| ). The first claim above is equivalent to: 



(B.32) 



/2£2 

-52 



(/2 - /o)52 

50 - 52 
1 - ho/h2 



and 



^ I - /lo//l5 
_/o - h 

\ hilib - go) 
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which, in turn, is equivalent to: 



J2 = 92 , 92 = /i2T 7-, (B.33) 

90 - 92 ho - /i2 



and 



^5 = h ^, /s = 95 —■ (B.34) 

/o - /5 50 - 55 

All these relations easily follow from the equations of the fgh-system. For instance, to check 
the first equation in ( |B.33| ), one has to consider the two elementary positively oriented triangles 
(3,3 + a;,3 + e) and + 1,} + e). Denoting the edge eu = (3 + w, 3 + e), we have: 

f2 + fi2 = fo + h (= -/i), fi292 = J290 (= h9i)- 

Eliminating /12 from these two equations, we end up with the desired one. This finishes the proof 
of ([Bl4|)-(pl^). 

Step 3. Checking the equations for the matrix D[i). Notice that the matrices 




L(e,0) = 

are upper triangular. We require that the matrices D{i) are also upper triangular: 

/ dii di2 di3 

D = i d22 d23 
V ^33 

It is immediately seen that the diagonal entries are constants. By multiplying the wave function 
^k,e,m{fJ') from the right by a constant (/i-dependent) matrix one can arrange that the matrices 
D(^) are traceless. Hence the diagonal part of D is parameterized by two arbitrary numbers. It 
will be convenient to choose this parametrization as 

(dii, d22, ^33) = ( - (2a + /3)/3, (a - f3)/3, (2/3 + a)/3) . 

Equating the entries 12 and 23 in (|B.17D -( ^J9|) , we find for an arbitrary positively oriented edge 
e = (31,32) G-E(T£): 

^12(32) - c?i2(3i) = id22 - du)f = a(^u{i2) - u{ii)^, 

^23(32) - ^23 (31) = (dss - d22)g = f3(^v{}2) - vinfj- 

Obviously, a solution (unique up to an additive constant) is given by 

di2 = au, d23 = Pv. 

Finally, equating in ( [B.l?| )-( ^.ig| ) the entries 13, we find: 

^^13(32) - c?i3(3i) = ^23(31)/ - (^12(32)5 (B.35) 
= Pv{n)(^u{i2)-u{u)^ -au{i2)(v{i2)-v{n)y (B.36) 
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Comparing this with ( [4.16| ), ( |4.18| ), we see that ( 5.22| ) is proved. 

Step 4. Equations relating the matrices Ck/^m and D(^). It remains to consider the 
equations ( B.2C| )- (|B.22 ). Denoting entries of the matrix C by Cij, we see that these matrix equations 
are equivalent to the following scalar ones: 





Cl2 + di2 = 0, 






(B.37) 




C23 + d23 = 0, 






(B.38) 




Cl3 + di3 = 0, 






(B.39) 


(C33)fc+l/,m 


— {Cu)k,e,m + ^33 — 


du 


= 1, 


(B.40) 


(C33)fc/+l,m 


— {cu)k,e,m + ^33 — 


du 


= 1, 


(B.41) 


(C33)fc/,m+l 


— (Cll)fc/,m + (i33 — 


du 


= 1. 


(B.42) 



(In the last three equations we took into account that du, ^33 are constants.) It is easy to see that 
the equations ( B.37D , ( |B.3S|) are nothing but the constraint equations (|5.5| ), (|5.6|) , respectively. We 
show now that the remaining equations ( B.3S| )-( |B.42 ) are not independent, but rather follow from 
the equations of the fgh-system and the constraints ( B.37| ), ( B.38 ). We start with the last three 



equations, and prove the claim for (B.40), since for other two everything is similar. As in the proof 
of Proposition we write the formulas here with m = 0. Writing ( B.40 ) in long hand, using the 
ansiitze (15^201 ), (|]2|), we see that it is equivalent to 

1/ho , l/ho 



a 



/5 



(k + l) 



fo9o + 90 f 3 + h93 

1^5 



k- 



fo9o + 90/3 + /sffs 



/252 +52/5 + /555 /252 + 92/5 + /555 



(B.43) 



But this follows immediately from (|B.6| ), (B.32). Finally, we turn to (B.39). Actually, since the 
entry 13 of the matrix D is defined only up to an additive constant, this equation is equivalent to 
the system of the following three ones: 



(Cl3)fc+1/ 

{Cl3)k,e+l,m 

{Cl3)k,e,m+1 



{Cl3)k,e,m 


+ di3 


- di3 


= 0, 


{Cl3)k,e,m 


+ di3 


- di3 


= 0, 


{cu)k,e,m 


+ di3 


- di3 


= 0. 



(B.44) 
(B.45) 
(B.46) 



As usual, we restrict ourselves to the first one. Upon using the equation (B.35) and the constraints 
( B.37 ), ( B.38 ), we see that it is equivalent to 



{C-I3)k+I,£,m {Cl3)k/,m +5o(ci2) 



k+l,e,m 



foic23)kJy 



0. 



(B.47) 



Writing in long hand, in the representation with m = 0, we see that the terms proportional to fc+ 1 
and k vanish identically, while the vanishing of the terms proportional to i is equivalent to: 



1 



1 -90/95 



1 



1 - /0//2 



/12/15 f 292 + 92/5 + f595 ^2/i5 f292 + 92h + fm 

But this follows immediately from ( [B.32| ) and the formulas 



55 



90-95 

hn - /15 



92 



90 - 92 

I 

■ /lO - /l2 



which are similar to (and follow from) the equations ( B.34 ), ( B.33 ). 
This finishes the proof of Theorem Il9|. I 
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